Minimal conductivity of rippled graphene with topological disorder 
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We study the transport properties of a neutral graphene sheet with curved regions induced or 
stabilized by topological defects. The proposed model gives rise to Dirac fermions in a random 
magnetic field and in the random space dependent Fermi velocity induced by the curvature. This 
last term leads to singular long range correlated disorder with special characteristics. The Drude 
minimal conductivity at zero energy is found to be inversely proportional to the density of topological 
, disorder, a signature of diffusive behavior. 
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I. INTRODUCTION 




Since its experimental realization [3, 0, S] , graphene has been a focus of intense research activity both theoretical 
and experimentally. The origin of this interest lies partially on the experimental capability of explore the transport 
properties that show a set of interesting features related to disorder and which open the way to graphene based 
r . electronics (4|. 

One of the most intriguing properties of graphene is the observation of a minimal conductivity at zero frequency 
in undoped suspended samples that in early measurements was argued to have a universal value of the order of e 2 /h 
independent of the disorder concentration and a factor of n bigger than the one predicted by theory 0, H, 
m m m q ei Il6j |. Other experiments in both mechanically deposited and graphene grown on a substrate 
1 have found a bigger dispersion in the coefficient of the universal behavior while more recent calculations 
13 . l2fj have casted some doubts on the disorder dependence of the numerical coefficient and the actual situation 
remains unclear. 

Another peculiarity of most of the graphene samples is the existence of mesoscopic cor rug ations (2ll . [22l . |23| whose 
possible influence on transport properties only now starts to be explored [2l|, [H, [H, [2(| [2?| ■ Although the observed 
l— "™ 1 ■ ripples were invoked from the very beginning to explain the absence of weak localization in the samples @, l28l |29| . 
I | there have been so far few attempts to model the corrugations based either in the curved space approach with [3Ct l3l| 
■ or without topological defects [32j , or on the theory of elasticity [U, [H, [H, [35[ . The possible physical implications 
of the ripples in connection with the charge anisotropies has been revisited in very recent works [3|| [37j • 
0^ ■ In refs. [3(1 [M| we proposed a model for rippled graphene based on the presence defective rings (pentagons and 
heptagons) in the samples. These types of topological defects have been observed in nanotubes and in bombarded 
graphite and are known to be a natural way to get rid of tensions in the hexagonal lattice [HJ . It is then natural to 
0^ " think that some of the defects that were either present in the graphite sample or formed during the very energetic 
| procedure of mechanical cleavage, stay quenched in the two dimensional samples. 

It is clear by now that the low energy electronic properties of graphene are very well described by the massless 
Dirac equation in two dimensions, a fact coming from the symmetries of the hexagonal lattice and also obtained in 
the tight binding approximation [391 ]. It is also known that the Dirac points are very robust to deformations of the 
lattice [4(| so we have modelled curved graphene assuming that the Dirac points are not affected by the presence of 
ripp les and hence that the curved samples can be described by writing the Dirac equation in the given curved surface 
[4lf . Within this formalism studied recently the electronic structure of the sample with topological defects [3(1 HH 
and that of graphene with smooth curved regions [HJ . In the last work we emphasized the fact that curvature gives 
rise not only to an effective magnetic field, a property known from the early times of graphene [42], [H, [44j, but also 
to an effective position dependent Fermi velocity which can have strong influence on the physical properties of the 
system. 

In this work we continue studying the physical properties of curved graphene within the "gravity" approach. The 
electronic properties were explored by means of the two point Green's function of the electron and we could make 
advances keeping a fixed number of defects at given positions of the lattice. To study the transport properties is a 
much more difficult task. We need to assume a density of defects with some statistical distribution and to average 
over defects. We apply the standard techniques of disordered electrons (45j | to rippled graphene by averaging over the 
random magnetic fields and over the effective Fermi velocity induced by curvature. We will make special emphasis on 
the case of having topological defects. The smooth curvature case can be implemented easily. We find that averaging- 
over the space-dependent Fermi velocity treated as a random scalar field in the case of having topological defects, 
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gives rise to a singular, long range correlated interaction that affects severely the one particle properties of the system. 
We compute semiclassically the zero frequency conductivity and find that it depends on the inverse of the density 
of disorder, a behavior characteristic of diffusive systems. We argue that even if the Drude value obtained in this 
work gets renormalized by quantum corrections to the universal minimal conductivity there will still be a region in 
parameter space where this model differs from the ones studied previously. 

The article is organized as follows: In Section |TT] we review the model of topological defects and establish the 
effective Hamiltonian. In section ITO1 we define the statistical properties of the fields induced by the defects and apply 
the replica trick to get the effective four Fermi interaction. We see that the interaction is anisotropic and singular in 
the forward direction, a behavior due to the long range character of the conical singularities. We discuss this issue 
and extract the interaction coefficients that will appear in the computation of the lifetime and the density of states 
of the disordered system. Section IIVI contains the semiclassical expansion of the sigma model to get the low energy 
behavior of the system. We discuss two possible saddle points and compute the effective potential to establish the 
symmetry breaking minimum, a basic point to the rest of the calculation. From there we deduce, after some lengthly 
calculations detailed in the Appendix A, the diffusive behavior of the system. Section IVl contains a discussion of the 
results and open questions. 



II. THE MODEL 

In this section we follow closely refs [3(J [H| to describe the explicit form of the potential generated by the defects 
and its statistical properties. The model for corrugated graphene is based on the presence of defective carbon rings 
at arbitrary positions in the lattice. It is known that substitution of an hexagon by an n-sided polygon with n greater 
(smaller) than six gives rise to locally curved portions in the sample with positive (negative) curvature. It was argued 
that the presence of an equal number of heptagons and pentagons would keep the sample flat in average and mimic the 
corrugations observed in the free standing samples. In the mentioned references, we studied the electronic properties 
of the sample by computing the electron Green's function in the curved space defined by a given number of heptagons 
and pentagons located at fixed positions of the lattice. In order to study transport properties we need to consider a 
density of defects located at random positions. 

The behavior of the electrons in curved graphene is described by the Hamiltonian: 

H = iv F J d 2 r^^ (0„ - fi^r)) V- (1) 

The curved Dirac matrices 7 M (r) are related with the usual constant matrices 7° by 

7 ^(r)= e ^(r) 7 Q , 

where is the tetrad constructed with the metric tensor. It is this factor in combination with the determinant of 
the metric y/g gives rise to a space dependent effective Fermi velocity. In terms of a tight binding language this 
term can be modelled as a global modulation of the nearest neighbor hopping [3(| or of the average distance between 
carbon atoms [37j induced b y c urvature. The term contains the spin connection and the possible extra gauge 
fields induced by the defects j38|. It is given by 

Q^(r)=r /J (r)-r a A«(r) (2) 

where T M is due to the spin connection and the non-abelian part t q A^ is related to the holonomy and will be discussed 
later. 

In (3^, [3l| we described the curved space generated by a arbitrary number of topological defects by the metric 

9ij = eA( % (3) 

where the conformal factor A(r) takes the form 



N 

\ " H 



a* 



where fXj is a constant related to the defect (or excess) angle of the disclinations and a* is a constant of the order of 
the lattice spacing, interpreted as the radius of the "core" of the defect. The specific form of the conformal factor Q 
gave rise to strongly diverging one particle properties like the local density of states. We will see that in the procedure 
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of averaging over disorder one particle properties will remain singular while two particle properties, like the Drude 
conductivity, will be convergent. 

It is well known that the presence of topological defects has other consequences besides of curving the graphene 
sheet. Odd membered rings can mix the Fermi points that can be implemented by a nonabelian gauge field [42. |43|. 
Also if various defects are present, an extra phase appears due to the non commutativity of the holonomy operators 
associated to the valley mixing phase and the proper Berry phase that fermions get when they surround the group of 
defects (46l \at\ . All these phases are naturally incorporated in this formalism by generic external nonabelian gauge 
fields A a (r) in equation fTJ) . This also accounts for the disordered graphene classification described in [14j, |48| . In the 
case of having an equal number of pentagon and heptagonal defects it can be shown that only the (abelian) gauge 
field associated to the conical singularity remains [381 ] we can restrict ourselves to the scattering problem around a 
single Fermi point. We will comment on the possible modifications of the calculation that a more general case would 
induce in section fVl 

The value of \fij\ = /i is 1/24 for both pentagon and heptagon rings. We use [i as a perturbative parameter around 
fiat space and expand the determinant of the metric g(r), the zweibeins e^'(r) and the spin connection il(r) in (fT]). 
To first order in fi the Hamiltonian is 

H = iv F J d 2 r[(l + A(r))^d^ + ^ 7 l (cU(r))# (5) 

Eqs. (J4j) and ((5j) are the basis of our calculation. We have a disorder-induced modification of the Fermi velocity 
described by A(r) and a gauge field- type term if we identify A4 ~ 9jA(r). Notice that although the gauge potential 
Ai seems to be a total derivative, the effective magnetic field is not zero due to the singular form of the function A(r) 
(U]). Also due to the presence of the conformal factor A(r) in the effective Fermi velocity of the Hamiltonian ([5]) even 
if the potential could be be gauged away by means of a singular gauge transformation (a coordinate choice in our 
case), its effects on the observable quantities will show up as coming from the other term. 



III. AVERAGING OVER DISORDER 



We will study the transport properties of this topologically disordered graphene by assuming a random distribution 
of an equal number of pentagons and heptagons. The statistical properties of these topological defects where analyzed 
in part in 49] . It was proposed there that the gaussian disorder induced by a random distribution of topological defects 
can be described by a single dimcnsionlcss quantity A proportional to the average fluctuations of the -non abelian- 
vector potential representing a vortex at the position of the defect. 

(A(r)A(r')) = A<5 2 (r-r') (6) 

The infrared behavior of A was analyzed in [4!| and shown to diverge logarithmically with the size R of the system 
in the case of unbounded disclinations 

A = 2^ log (^) (7) 

while it remains a of constant value 

A oc $ 2 n disl b 2 (8) 

in the case of having pentagon-heptagon pairs bounded into dislocations with density ridtsi and with an average 
distance b. 

We will work in the intermediate case of having a random distribution of an equal number of five and seven rings 
at average distances d such that the "volume" occupied by virtual strings pairing the defects is small compared with 
the total size of the sample (d/L) 2 <C 1. 

The most important issue in this work is the new random field A(r) associated to the Fermi velocity modification 
given in eq. (|3]). We will assume for this scalar field a zero mean value (A(r)) = and 



(A(r)A(r')) = na* 2 log 



r — r 



(9) 



where n is proportional to the areal density of defects and a* is of the order of the lattice spacing. This correlator 
diverges in both the infrared and the ultraviolet limits. It induces over the random magnetic field Aj(r) the average 

{A i (r)M*')) ~ na* 2 5 t3 5{v - r'). (10) 



4 



The origin of the correlator ([9]) can be understood when we consider the nature of the defects. In a geometrical 
description of defects in two dimensional crystals 50] , the equations of motion for the metric tensor (r) reduce to 
a unique equation for the conformal factor J4}: 



A' 



V 2 A(r)=^g-«5(r-r i ). (11) 

3 

We can rewrite ([J} in terms of the Green's function of pip : 

N 

A ( r ) = E H J dr ' S ( r ' r ^( r r ')- (12) 



where the Green's function has de form 



for r ^> a* , what justifies equation ©. 
In momentum space the correlators are 



K(t - r') « log 



(13) 



(A(p)A(-p)) = — (14) 

(A i (p)A j (-p))=na* 2 6 ij (15) 

Notice that the A term describes a new type of disorder: Dirac Fermions in two space dimensions with a random 
velocity, a problem that, to our knowledge has not been addressed in the early literature [5ll l54j . 

In momentum space representation the action corresponding to eq. JSJ) reads: 

S = \J dk^j (luj° — vpjftj ip — —vp J dk—ipjA(k)^j ^- J dpdkA(p)ipjkip, (16) 

where dk = and 

A(P)= / f^log|-|. (17) 
J Zir a 

We integrate out A using eq. I]14p as a quadratic action. It is interesting to note that this term is similar to the 
interaction between curvatures in a continuous hexatic membrane where l/na* 2 plays the role of the hexatic stiffness 
constant (Hoj . 

Replicating the fields and integrating out A in (|16[) we get: 



S = So + v 2 F ^J dk dk'T{k, fcOWv7Va)(^7*^), (18) 

where summation over replica indices a, b is assumed, A = 27r/i 2 na 2 is a dimensionless parameter proportional to the 
density of defects n and 

The constant term in the interaction vertex (|19p comes from the random magnetic field and would give rise to the 
standard result found in previous works. The term coming from the Fermi velocity is anisotropic and singular when 
k — * k', a signature of the infrared singularity associated to the effect of a conical defect at infinite distances from the 
apex. Because we are dealing with elastic scattering, the modulus of the momentum is conserved provided that the 
energy uj in the process is conserved. The function T(k, k') in (|19[) can be written as a function of the difference of 
the scattering angles, <f> = 0}. — 6y. 



cos 2 (</>/2) 1 
sin 2 (0/2) ~ 4 
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This function diverges at scattering angles (f> = (forward scattering) and the one particle properties of the system 
like the density of states or the imaginary part of the self-energy will show anomalous behavior when compared with 
their counterparts in short ranged scattering processes. The divergence in ([20| can be regularized by including a 
cutoff 6: 



.Mn 2 (</)/2) + 6 2 

The meaning of the cutoff can be understood from the origin of the correlator for the function A(r) in eq. ©■ Instead 
of having a long rangi 
correlator of the type 



of having a long ranged propagator corresponding to infinite range defects that behave like p., we may consider a 



what corresponds in real space to changing (| 13[) by a modified Bessel function of the second kind: 

K(r-r') =K (<5|r-r'|). (23) 

The leading term in the expansion of (f2"3"|) for small 6 is 

K(r-r') ralog(£|r-r'|). (24) 

Using eq. (f2"2"|) instead of eq. (|14[) we arrive at eq. (|2"Tj) , where we have redefined <5 as 6 — 8/k 2 . We can assume that 
the momenta involved in the problem are of the order of the wavelength A of the states near the Dirac points (or the 
localization length if those states are localized) so we can consider that 6 is of the order of X/x where x is the biggest 
length scale of the problem comparable to the system size L. We will make the important assumption that x > ^- 
We will see that the Drude conductivity is independent of 6 and then well defined. 
The free action in (fT8|) is the usual: 

So = / dkiPaiuj-f - v Fl k + iriM)ip a , (25) 



The replica indices run over 1 to 2N . The first N indices are associated to advanced fields (+irj in the free action), 
and the second ones to retarded fields, related to —it], with the obvious definition of the matrix M: 



( ° ) ■ <-><>> 

' 2N 



This ?7-term breaks explicitly the 0(2N) symmetry of the action (fTH)) leading to possible massless excitations [55|. 

In addition to the forward divergence the scattering function ([21]) is highly anisotropic [56| what will lead us to 
deal with a multichannel scattering problem in contrast with the sort ranged scattering case where the scattering only 
occurs in the s-channel. We will see that despite the anisotropic scattering the diffusion process will have an isotropic 
behavior described by a scalar diffusion constant. 

In order to keep track of the role of each channel, we decompose T(<p) in harmonics: 

r(^) = ^r„ e -"^. (27) 

n 

We will see later that only the n = 0,±1 channels will play a role (notice that T n = r_„). Their explicit values as a 
function of the cutoff 5 are: 

r = I - 5/4 + 0(6) , T ± = 1-2 + 0(6). (28) 



IV. THE SADDLE POINT APPROXIMATION IN THE NONLINEAR a MODEL AND RESULTS 

In this section we will construct the low-energy field theory describing the large scale behavior of the system (i.e., 
at length scales larger than the elastic mean free path, I) following standard procedures [45j]. 
The interaction term in (|18p can be written as 
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Sim = \^IY1 I (dk)(dk')r nX nCk)xUk')&af^a)(j b ^ b ), (29) 



where Xn(k) — e m0h . Now we proceed to simplify the model (|29p . The current-current interaction in ([29")) can be 
transformed into a density-density type term [571 ] : 

Sint = ~ Kff^l ^ / dkdk'r n Xn(k)xl(k')[(^a^b)(^b^a) ~ (^al5^b){^bl5^a)}, (30) 



where we have redefined the coupling constant as g = — 2A. 

The four fermi interaction in (|30[) can be decoupled by means of a Hubbard-Stratonovitch trnasformation 

S = S + f dk^-^i^n {Ql+WD+i^nXnik^iQn-l!^)^ 



(31) 



where Q n and II„ are Hubbard-Stratonovitch fields. Because we have defined retarded and advanced fermionic fields, 
Q and II will be 2N dimensional matrix fields. 

We can further simplify the calculations by performing a unitary transformation that diagonalizes the fermionic part 
of the action in (f5Tj) . and leaves unchanged the functional integration measure: ip — > Uip, $ — > i^U + , Q n — » UQ n U + , 
n„ — > C/II n {/ + , 75 — > U^i§U + . Without loss of generality we will name those transformed fields as the original ones. 
Integrating out the fermionic modes leads to the usual final form for the non linear a model: 



S eff = [ TrlnG- 1 -L- £ r„(Q* + l£), 



(32) 



where the Green's function is defined by 

G- 1 = Gq 1 + irjM + i T nX n(k)(Q n - h 5 U n ), (33) 



and 



i i w + vpk \ 



w — upfc 

Next we will make a saddle point approximation and seek for a solution of the saddle point equation 

S Seff 



5{S n ) 



0, (35) 



In contrast with the case of isotropic short ranged scatterers 0, where a single equation is obtained, ([35]) represents 
an infinite number of coupled equations: 

^ ^ = f W ^t-i -^ n l k) /fwov - ( 36 ) 
2 3< V G 1 + i r„ X n (fc) (5„) 

Following [59], in the limit w — > 0, and despite of the fact that the scattering mechanism is anisotropic, we can try to 
find a solution of the type 

{Qn) = fS n0 M. (37) 

Of course, the solution of the mean field equations is not unique, and other vacua with different properties may be 
found. The value of the trial function / is: 

/=^(e^-lp, (38) 
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where K is an ultraviolet cutoff. Now we are free to make a choice over the mean field values of Q and EL The 
standard choice is (Qo) = (So) an d (EIq) = 0. With this solution for (Q) we arrive to the one particle relaxation time 
or lifetime 

L= VF K(e^ (39) 

which turns out to be a constant whose dependence on the disorder strength is typically non-perturbative: 
r -1 (g) ~ exp(l/g) . This value for the lifetime depends on both the ultraviolet cutoff K, and on the infrared cutoff 5 
through the scattering coefficient Tq ^ 4 — 5/4. The ultraviolet cutoff W can be removed by renormalization group 
techniques as done in [58| but the infrared cutoff S remains and the one-particle properties of the theory will depend 
explicitly on it 1 60] . We will see later that the transport properties ( i.e. two-particle properties) are independent of 
such regulator [681 ] . 

We can now compute the averaged density of states at the Fermi level, 

p(0) = -^Im/g<7«(*), (40) 

the index N in the denominator is the replica index and the limit N — > will be taken at the end of the calculation. 
From now on, we drop any reference to this limit, bearing in mind that it must be taken when computing observable 
quantities. 

The density of states as a function of the lifetime (f3"9"| is: 

p(0) = 4rf- (41) 

Before proceeding to compute the quantum fluctuations around the chiral symmetry breaking solution of the saddle 
point equations we will make a comment on the solution (S n ) = 0. This solution leads to a vanishing density of states 
at the Fermi energy and to a sublinear frequency dependence: p(to) ~ |w| Q , with a being a function of the strength 
of the disorder, a behavior reported in (5l| . In order to determine the true minimum of the Q field action we have 
computed the effective potential as a function of / defined as in ([55)1 . This effective potential is the same as that of 
the Nambu-Jona Lasinio model [6l|, [62|, [63| being a function of / defined as in ([55)) , when (Hq) is taken to be zero: 



The result is well known. Minimizing this potential the two possible solutions are again the solution / = corre- 
sponding to a zero value for the effective potential and the symmetry breaking solution / = Ke x l 2 ~ K ^ 2 e^ giving 
V e ff = — -f^-/ 2 - This broken symmetry solution equivalent to (|38[) is then a minimum of the theory. 

We will then proceed computing the physical properties of the quantum field model built around the broken 
symmetry solution. We will see that the physics obtained for this case is typically non-perturbative irrespective of 
the strength of the disorder. 

The technical details of the rest of the computation are given in appendix A. Once we have calculated the value of 
the leading configuration of Q n from the saddle point equations (|35[) . we expand the action Q32p around this value, 
setting Q n = (Q n ) + 5Q n and retain in the expansion terms up to second order in 5Q n : 

S*{S) + 5Q n * S * 5Q m + ... (43) 

The * means that the functional derivative is evaluated in the saddle point solution for Q n and II„. 
The ultimate goal is to compute the action for the massless modes: 



SS = / dqSQo^-^- (r) + Dq 2 ) SQ , (44) 



(2r)r , 2 
- — — [r] + Dq 

from where we can extract the diffusion constant is D. From eq. (|A18|) we get 
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where the coefficients To and Ti are given in (|28[) and their difference is Tq — T± — § + 0(5), hence the diffusion 
coefficient is well defined when the cutoff 5 is send to zero. 

Finally we can compute the semiclassical value for the static DC conductivity using the Einstein relation for the 
two diffusive channels SQo and SHo in eq. ([50)1 : 

e 2 4e 2 4 

aDC=4 - p (0)D = -2— (46) 

The factor of 4 comes from the spin and valley degeneracy and the factor 2 comes from the two diffusive channels. 

The result in (|46[) depends on the coupling parameter g — Aw^na 2 , which contains information about the type of 
disorder and the density of disorder n. 



V. DISCUSSION AND OPEN QUESTIONS 



In this work we have addressed the effects of curvature on the transport properties of corrugated graphene sheets. 
We have shown that coupling the Dirac field to a curved space gives rise to an effective potential whose general form 
and statistical properties depend on the metric. The main feature of the geometrical description is the appearance 
of an effective space dependent Fermi velocity which gives rise to a random scalar field coupled to the kinetic energy 
term in the Hamiltonian. 

Smooth curved regions in graphene give rise to standard short range correlated disorder as the one studied in the 
literature [?J • The presence of topological defects in the sample either as the main source of curvature or as a way to 
stabilize the ripples in the mechanically deposited samples gives rise to singular, long range correlated disorder that 
leads the system to an intrinsically diffusive behavior even at the neutrality point and without the need for external 
scattering sources. 

We have studied the behavior of the conductivity at zero voltage of a graphene sample in presence of a finite 
density of topological defects. The result given in eq. (|46[) shows that the conductivity is proportional to the inverse 
of the density of defects, a distinctive feature of a diffusive system. A similar non universal behavior has been found 
in the same system and has been attributed to the effects of random coulomb scatterers present in the substrate 
[l9l [20l . [gil, [H| . A crucial difference is that in the mentioned references the graphene sample is either heavily doped 
or it has a nonzero carrier density due to a local field effect induced by the Coulomb impurities. In our work the 
density of states is generated by the disorder as in refs. (7I [Toj. 

Another noticeable feature of the model presented in this work is the strong dependence of the one particle properties 
on the parameter 5 regulating the infrared behavior of the model. The situation here is even worse than that of a 
two dimensional electron gas in a long range correlated random magnetic field discussed in ref. [59l ]. There the one 
particle relaxation time was found to diverge as the infrared cutoff is sent to zero but the finite density of states made 
the transport relaxation time r tr finite. In our case Tt r also depends on the density of states at the Fermi level but 
now the DOS is divergent for <5 — > 0. As we have seen despite this singular behavior we have obtained a finite Drude 
conductivity due to the particular dependence of the diffusive constant with the density of states (eq. I|45p). 

In section III we have introduced the parameter S defined in a phenomenological fashion as the ratio between the 
characteristic length scale of the defect x, and the wave or localization length A of the states around the Fermi energy. 
In a semiclassical approximation to the problem this parameter is essentially uncontrollable. We can nevertheless made 
an estimation of the range of applicability of our results by assuming that the localization length can be obtained 
from an analysis of the quantum corrections to the conductivity. The diffusive regime is characterized by a static 
mean free path I = VFT tr greater than the localization or wavelength A but smaller than the system's size. The mean 
free path can be estimated using expressions (|59"|) , (|4"Tj) . and (|43|) and assuming that a* ~ a and \ ~ L. We thus find 

j. 2 m V2 1 



3// 



(W^ 



77 

imp 



from where we can get a lower bound for the density of defects in the case A < I : 

4 L 
V A 5 

In the same spirit, an upper bound can be estimated using the condition I < L: 

4 I 
9/? ~L\ 



n lmp < ( 48 ) 
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In refs. [13J, [66| the possible fixed points of the total conductivity where classified according to the symmetries of the 
original Hamiltonian in a renormalization group analysis. As our disorder term preserves both chiral and time reversal 
symmetries, the final conductivity once quantum corrections are taken into account should flow to the universal value 
of 4e 2 /7r/i. We note that in previous works this universal value is already obtained at the Drudc level. The topological 
disorder discussed in this work sets as an initial condition of the RG flow a rather different - disorder dependent - 
value that can - or not - flow to the usual fixed point. The analysis of the quantum corrections to the conductivity 
(HI]) is beyond the scope of this work and will be reported somewhere else. 
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APPENDIX A: SIGMA MODEL CALCULATIONS 



Once we have calculated the value of the leading configuration of Q n from the saddle point equations (|35|) . we 
expand the action (|3"2")) around this value, setting Q n = (Q n ) + SQ n and retain in the expansion terms up to second 
order in SQ n : 

Sn{S)+5Q n * S * SQ m + ... (Al) 

oQnOQm 

The finite density of states of eq. (|41 [) allows us to use the usual identity in the integration of momentum derived for 
the two dimensional electron gas: 



J — * p(0) J ' de k , v F k 



As a consistency test, computing the density of states with this change we get the expression (l4Tj) . This change of 
variables will simplify the calculations. The * in (|A1|) means that the functional derivative is evaluated at the saddle 
point solution for Q n and II n . Also, other terms of quadratic order appear in (|Alj) with crossing functional derivatives 
in the fields 5Q and 511. We will see shortly that these derivatives are zero and the former fields are not coupled. 
These derivatives are 

<5 2 S* I 

OQnOQm 

and 



l - FnSQl + \J (dp)(dq) Y ^nT mXp (k)Xm( P + q)SQ n TrG(k)G(p + q)SQ m , (A2) 



SUn Sn n Sn m m ' n = ~l^ r " m " ~\j {dp){dq) 2 T ^mXn{p)Xm{p + q)6U n Tr l5 G(p) l5 G(p + q)6U m . (A3) 



Since the spectral functions are peaked at the Fermi energy, we can effectively restrict the q integration to values 
around the Fermi point, q << Kp and write Xm(p + q) ~ Xmip) causing an error of the order of 0(q). In (|A2[) the 
integral will be dominated by the product G R G A , i.e., by the off diagonal sector of the fields 5Q^ leading to the 
diffusive pole behavior for the fields Q n . By contrast, in (|A3j) using the symmetry property 756*^75 = — G A we can 
see that the diffusive pole will come from the product G A G A and its retarded-retarded counterpart 5Il ++ ~~ will 
be the channel for the diffusive behavior (also, the minus sign appearing in this symmetry corrects the relative sign 
between the second terms in (|A2j) and (IA3jl ). 

Let us define the quantity C mn (r), q) as 

C mn (v, l) = \j (d P )Xn(p)x m (p)TrG R (p)G A (p + q). (A4) 

To obtain an action for the modes Q n and n„ at small q and 77, we expand G R up to first order in 77 and to second 
order in q in equation (|A4|) . 
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The first term in this expansion , without any reference to the replica index, is 



C° nm = \J (dp)e^ n+ ^ e Tr ( 



\(v&B 
^-L-^ J I 1 



We have denoted e = 77 4- In terms of these variables, we have 



4 ^+27 

Expanding (|A5j) up to first order in 77, and using l|4ip. we get 



C° nm = ~p(0)6- nm — Lj-. (A5) 



= ~ \p(0)5- nm (2r) 2 Ti + 0(, ? 2 ). (A6) 



The constant term in (|A6j) coincides with the term proportional to 8Q\ and <5II 2 in (|A2j) and (|A3j) respectively. This 
mass contribution to the action is 

^I^tE ^ - l) (*£ + ■ ( A7 ) 

We immediately see that only the modes SQo and <5IIo are massless, and they will responsible for the diffusive behavior 
of the system exactly as happens in the 2DEG 59] . In what follows we will eliminate the II field. It represents another 
diffusion channel that does not mix with the Q's and plays the same role. We will simply multiply by two the final 
result. The next term in the q expansion reads 

C« m = I / (dp)e i ^ B v F q cos 6 (- .7/ — V2 + -, • (A8) 

4 J \ [vfP — te)(VFP + it) (—VFP~te)(—VFP + ie) J 

Note that in the case of short ranged isotropic scattering this term vanishes due to the angular integration. In our 
case, however, the presence of e l i m + n ) e allows linear terms in q, coupling the massive modes 5Q±\ to SQo- The integral 
in (|A8|) is performed changing to the energy variable and noticing that the angular integration gives a non-zero result 
only when n = —m ± 1: 

1 f 1 

Cl m = -v F q(5nm~i+S nm+1 )p(0) / de p - — _ , , , (A9) 



or, after setting 77 = 0, 



CL = -^P(0)(<5„ m -i + <Wi)(2r) 2 . (A10) 



or, rearranging signs, 



To compute the next term in the q expansion, we will use the following trick. The trace of the product of Green 
functions in (|A4p can be explicitly written as 

TrG R G A = ( VF P~ le i \ ( »Hp+q|+« Q | ^ ( A11 ) 

\ —VFp—ie/\ — uplp+ql+^e/ 

TrG R G A = — — —. + — -. (A12) 

vpp — it vf |p + q| + it vfP + it vf |p + q| — it 

We immediately see that the second term in the right hand side is the complex conjugate of the first term, thus, 

TrG R G A = 2Re ( : —. ^ . (A13) 

\VfP - it v F |p + q| + it J 

The terms already calculated in the expansion of C nm (7?, q) can be easily derived with this trick, but where we make 
a real profit of this simplification is in the calculation of the term q 2 : 

Cl m = \fRe I (^c^-J— ( r ^- 3 - - . (A14) 

2 J vfP" it \ (VFP + it) 2vfP(vfP + it) / 
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If we compare the angular part of (|A14p with the corresponding part in (|A8[) we see that after performing the integral 
in (|A14|) there are terms of the type <5_ nm ±2 together with terms 5^ nm which generate couplings between the zero 
modes SQo and the massive SQ±2, and SQoSQo respectively. The couplings involving 5Q±2 being of order q 2 will 
produce terms of order q A and can be neglected. We will only keep the terms independent of 8 in (|A14[) . from which 
we will extract the diffusion coefficient D for the massless diffusive mode SQq. The result for (|A14[) only taking into 
account the terms proportional to b- nm is (we shift the pole at vfP = in the second term in the integrand and take 
the real part, the first term will not contribute to this real part): 

C 2 nm = 4^-™p(0)(2t) 3 . (A15) 
Collecting all the terms, the action for the modes 8Qo and SQ±i is: 



J F „±1 

Y x T Q p{Q)v F q{2 



,2 



t (<5Qo<5Qi + <WQ-i). (A16) 

In order to get a theory for the n = modes, we integrate out the n = ±1 modes in (|A16|) . Using again (j4"Tj) we get 

, s ■ j d /J^im iQo ^gn + ^. (A17) 

If we simplify and use again (|4"Tj) we arrive to the final action for the massless modes: 

SS = f dq6Q ^P^ (rj + Dq 2 ) SQ , (A18) 
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